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Abstract 



Qualitative approach to homogeneous anisotropic Bianchi class A models in terms 
of dynamical systems reveals a hierarchy of invariant manifolds. By calculating the 

■ Kovalevski Exponents according to Adler - van Moerbecke method we discuss how 
\^ \ algebraic integrability property is distributed in this class of models. In particular 

■ we find that algebraic nonintegrability of vacuum Bianchi VIIq model is inherited by 
\ more general Bianchi VIII and Bianchi IX vacuum types. Matter terms (cosmological 

constant, dust and radiation) in the Einstein equations typically generate irrational 
^— ^ \ or complex Kovalevski exponents in class A homogeneous models thus introducing an 

■ element of nonintegrability even though the respective vacuum models are integrable. 

1 Introduction 

Einstein's theory of General Relativity provides a very elegant geometric picture of grav- 
^ ■ itational interaction. However, such a formulation becomes possible at the expense of 

making the field equations strongly nonlinear |||]. Even in rather symmetric case of spa- 
tially homogeneous space-times the dynamics of their general types (mixmaster models) 
turns out to be chaotic Q. 

On the other hand, growing importance of non-linearity in fundamental as well as in 
phenomenological physical laws and richness of structure induced by non-linear phenomena 
motivated strong interest in techniques allowing one to detect and properly describe the 
complexity inherent to non-linear systems. One such approach aimed at investigating the 
integrability of non-linear systems has been initiated in works of Kovalevski |^] who in 
her famous contribution to the problem of rotating rigid body considered the behaviour 
of the solution near the essential singularities in the complex plane. Similar intuition that 
the position of singular points (essential singularities ) depends on (arbitrary) integration 
constants and consequently that their mobility precludes the construction of first integrals 
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was formulated and developed by Painleve Ig. For the review of singularity analysis, 
Painleve property and integrability of nonlinear systems see [||. 

Yoshida developed precise methods of investigating the restricted case of algebraic 
integrability of homogeneous dynamical systems. In particular he obtained a necessary 
condition for algebraic integrability in terms of Kovalevski exponents (i.e. the nonexistence 
of irrational or complex KEs). Kovalevski Exponents § ^ allow for reconstructing a set of 
essential singularities for nonintegrable systems without need of constructing the solution. 
The usefulness of calculating the Kovalevski exponents has also been emphasized by Adler 
and van Moerbecke [^. They gave the criteria for algebraic complete integrability for 
Euclidean Toda-like systems. Pavlov Q essentially extended Adler-van Moerbecke formula 
into the form suitable for studying indefinite Toda lattices. Thus from the existence 
of irrational or complex KEs we shall deduce that the system does not have additional 
algebraic first integrals (i.e. besides the energy integral). Of course the nonexistence of 
such first integrals does not imply nonintegrability (see e.g. ||5|). On the other hand when 
the KEs are integer numbers we not only know that corresponding first integral exists but 
we are able to construct its form. 

In this paper we calculate the Kovalevski Exponents (KEs) in the wide class of ho- 
mogeneous Bianchi spacetimes represented as hamiltonian systems with indefinite kinetic 
energy form and exponential potential term. Our technique is a generalization of the 
methods previously applied in studies of Euclidean Toda systems ||8| . In the present paper 
we applied this method for the wide class of Bianchi models with matter and cosmological 
constant. Our investigations are complementary to those of Demaret and Sheen Q who 
performed the Painleve perturbative test for the Bianchi Type IX model with perfect fluid, 
radiation and cosmological constant in order to predict some probable chaotic regimes. In 
all cases but the cosmological constant case, the models studied do not pass the Painleve 
test, exhibit infinitely many movable logarithmic singularities and are therefore probably 



chaotic. In [10| (see also [|lT|) it was shown that the perturbation of an exact solution 
exhibits a movable transcendental essential singularity thus proving the nonintegrability. 
Further numerical evidence of the nonintegrability of the Mixmaster model by demonstrat- 
ing the existence of complicated dense essential singularity patterns and infinitely-sheeted 
solutions with sensitive dependence on initial conditions can be found in the paper by 
Bountis and Drossos [^] . 

It is known that there exists a hierarchy in the class of homogeneous cosmological 
models dictated by the dimensionality of isometry groups acting on constant time hy- 
persurfaces |15, 16]. This was in fact the strongest argument supporting the claim that 
mixmaster models represent a generic dynamical regime near the cosmological singularity 
1^, 16 1 . We shall explore in this paper how this structure is reflected in the property of 
algebraic integrability in the class of homogeneous anisotropic cosmological models. 



2 Kovalevski Exponents in brief 

In this section we shall recall the main facts concerning the Kovalevski Exponents. For 
this purpose, let us consider a system of flrst order ordinary differential equations: 



(2.1) 
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where the right hand sides Fi are rational functions of xi, . . . , Moreover, assume that 
the functions Fi are weighted homogeneous with weights Mi i.e. 

Fiia^'xi, . . .,a^"Xn) = a^^'Fi{xi, . . . ,2;„), 

where Mj = gi + 1. Then the system ( |2.1D is similarity invariant and there exists a 
particular solution 

where the constants Cj satisfy algebraic equations: 

Fi(ci, . . . ,c„) = -giQ, i = l,...,n. (2.3) 

When we consider variational equations of the initial system (^]^) about the reference 
solution ( |2.2| ) and substitute Xi = {ci + Ci)^~^* obtain the linearized system 

where 

The eigenvalues pi,...,p„ of the Kij matrix are called the Kovalevski Exponents. In 
iQ] it has been proven that a necessary condition for a similarity invariant system ( |2.l[) 
to be algebraically integrable is that all Kovalevski Exponents be rational numbers. In 
other words, if at least one Kovalevski Exponent is irrational or imaginary then the sys- 
tem ( |2.1| ) does not have algebraic first integrals. This result is central to our subsequent 
investigation. 

3 Bianchi class A models as indefinite Toda lattices 

Einstein field equations for homogeneous Bianchi class A models as well as for diagonal 
class B models can be cast into a hamiltonian form [^] 

W = ^ i-pl+Pp. +exp(4a)T/(/3+,/5_), (3.1) 



where 



y(/3+, /3_) = nl exp (-8/3+) + nj exp (4/3+ + 4\/3/3_) + 

+ nl exp (4/3+ - 4\/3/3_) - 2nin2 exp (4/3+) - 

- 2nin3 exp {-2(3+ + 2\/3/3_) - 2n2n3 exp (-2/3+ - 2\/3/3_). 

The variables a, /3+, /3_ have the meaning of expansion factor and anisotropy parameters 
respectively Q whereas ni, n2, parametrize the structure constants of isometry groups 
(defining the type of the spacetime). 

Hamiltonian (3.1) is a generalized indefinite Toda lattice i.e. it is of the type 

1 ^ 

'H = -{p,p) + ^Ci(iyi^{ai,q), (3.2) 

i=l 

where (• , • ) is a Lorenzian scalar product and (• , • ) is an Euclidean scalar product, a is a 
real vector. 
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4 Kovalevski Exponents for indefinite Toda Lattices 

Given generalized Toda lattice let us define new variables 

Ui = {ai,p), Vi = exp{ai,q). 

Transformation between {a,p) and {u,v) variables is a generalization of Flaska transfor- 
mation. 

Then one obtains an autonomous dynamical system with polynomial right hand sides 

N 

Vi = UiVi, iii = '^MijVj, Mij = -Cj{ai,aj). (4.1) 

The system (^]^) is of the type ( |2.lD and one can easily check that it is similarity invariant 
and possesses a particular solution (see (^^)): 

-El -Yl 

where the constants C/j, Vi fulfill the following algebraic equations (see (^.3|)): 

n 

-2Vi = UiVi, -Ui = MijVj. (4.2) 

i=i 

In order to discuss the single-valuedness of solutions one investigates behavior of so- 
lutions in the vicinity of particular solutions. Let us hence consider the first variation 

—Sui = 'S^Mij5vj, -rSvi = — — - -\ — ^-^r^. (4.3) 
dt ^ ^ dt t ^ ^ 

We seek for the particular solutions of ([4.3| ) in the form 

6ui = Citf-\ 6vi = r]itP-^, 

where ^j, r]i are constants and p are called Kovalevski Exponents and satisfy the following 
system of algebraic equations: 

n 

(p - 2 - Ui)7ii = Viii, {p - = ^ Mi.rij . (4.4) 

i=i 

The general procedure for extracting the Kovalevski Exponents is to solve the systems 
(P) and (0). 

Some specific cases can clearly be distinguished: 
{i)3i:V,^0;yi^j:Vj = 0, 

(ii) 3 ■.Vi.^O,Vj^O;y k^i,k^j ■.Vk = 0, 

(iii) 3 {i,j,k) : V ^ 0,Vj ^ 0,Vk ^ 0;y I ^ i,l ^ j,l ^ k : Vi = 0, 

and so on until the dimension of the system is reached. Integer Kovalevski Exponents 
are necessary for solution to be meromorphic. 
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In the case (i) pi = —1 and p2 = 2 are always solutions, other Kovalevski Exponents 
can be calculated from the generalized Adler - van Moerbecke formula: 



P 



3] 



where i ^ j, Mjj ^ and j is such that Vj ^ 0. 
In the case (ii) one should solve the system: 



P(P-I) 



/'(P-I) 



M,,, 



MijMji 



0. 



For this purpose it is useful to build an auxiliary matrix 



Ma M^j 



(4.5) 



(4.6) 



5 Vacuum Bianchi class A models 

In this class of model we have: 

al = (4,-8,0), a2 

al = (4,4,0), a5 = 

and 



4,4,4^/3 
4,-2,2\/3 



as 
ah -- 



4,4,-4^/3 
4,-2,-2^/3 



C = l^ni, 712,713, 

The respective Cartan matrix reads: 



-2?li7T,2, -2712723, -272i?l3) 



M, 



-48 





-nl 


-nl 


2711712 








—nl 


nl 


-nl 








2ni7i3 


— nl 


-nl 


nj 





2712713 





-nl 













711713 








-nl 


711712 





711713 





-nl 





71i?l2 


?^2n3 






(5.1) 



Application of the generalized Adler - van Moerbecke formula (whenever there is only 
one Vi 7^ 0) shows that all Kovalevski Exponents (corresponding to this case) are integer 
and equal: 

Pi = -1, P2 = 2, P3 = 4. 

Let us now check if Mij is degenerate. For this purpose we build an auxiliary traceless 
symmetric matrix, 

96711712 





-48nin2?i3 
-48?ifn27i3 



[det Mi- 







96ni7i2 


















96712713 







96712 Tig 

-487iin|n3 


4871i7l27l3 



96ni7l27l3 

An inspection of (^]^) allows us to formulate the following 
Conclusions 





96711712713 


-4871^712713 
-4871in27l| 





(5.2) 
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• for Bianchi I and Bianchi II models every block matrix is degenerate; 

• for Bianchi VIq and Bianchi VIIq models there is only one pair (Vi,V4) leading to 
nondegenerate matrix M14; respective Kovalevski Exponents read: B(VIo) pi,2 = 
i(l ± p3 = -1, p4 = 2 and for B(VIIo): z = p{p-l) = l± iV7; 

• from mixmaster models (i.e. Bianchi IX and VIII) there is six pairs of (V^, V^) leading 
to nondegenerate Mj^; respective Kovalevski exponents are the same for Bianchi IX 
and VIII models and read: z = 1 i: i\fl for pair (1,4) and z = ±2 for pairs (3,5); 
(4, 5); (4, 6); (5; 6). Because z = p{p — 1) then z = 2 corresponds to pi = —1, P2 = 2 
and z = —2 corresponds to pi^2 = ^(1 i i\/7). 

Corollary 5.1. Bianchi type II model is algebraically integrable. 

Corollary 5.2. Mixmaster models inherit their algebraic nonintegrahility from Bianchi 
VIIq model. 

6 Homogeneous models with cosmological constant, 
dust and radiative matter 

One can now consider cosmological models with nonvanishing matter terms. They are 
described by the hamiltonian: 

^ = \ {-pI+p}+ +Ph)+ exp (4a)y(/3+,/3_) + 
+ Aexp (6a) + pdust exp (3a) + Trad exp (2a), 

where A as usually denotes the cosmological constant and the subsequent terms denote 
other matter sources i.e. the amount of energy-momentum carried by dust and radiative 
matter, respectively. 



6.1 Models with cosmological constant 

Let us consider the case A / 0, Trad = /^dust = 0. In addition to six a, vectors in the 
decomposition (3.2) there exists now one more vector: 

«7 = (6,0,0) 

and subsequent formulae ( |5.1| ) and (5.2) are modified correspondingly. From the formula 
(4.5) one can calculate the following additional Kovalevski Exponents: 

Pi = 2 -2—— = - for 1 = 1,..., 6, 

(07,07) 3 

= 2-2^^4^ =3 for i = 1,2,3. 

{ai,ai) 

Other Kovalevski Exponents obtained from ([4. 61) are: 



Recall that z = p{p—l). Hence we obtained irrational and complex values of the Kovalevski 
Exponents. Their existence does not depend on the value of cosmological constant. This 
proves that the presence of A introduces an element of nonintegrahility in the model. 
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6.2 Models with dust matter 

Now we choose /Ujust 7^ and fix A = Trad = 0. Additional seventh Oj vector in the 



decomposition (3.2) is: 

a7 = (3,0,0). 

Adler - van Moerbecke formula (|4.5|) gives the following Kovalevski Exponents: 



Pi = 2-2— — = - for i = l,...,6, 

(«7,«7) 3 

/3j = 2-2— ^- = - for z = l,2,3. 

(aj,ai) 2 



Further Kovalevski Exponents obtained from ( ^.61) read: 

1 , \/59T 

Zi 2 = ± . 

20 60 

Therefore irrational and complex values of the Kovalevski Exponents exist irrespectively 
of concrete value of /j, or Bianchi type. 

6.3 Models with radiative matter 

We take Fj-ad 7^ and A = /^dust = 0. Additional vector in the decomposition ( |3.2| ) reads 
now: 

a7 = (2,0,0). 
Adler - van Moerbecke formula (|4.5|) yields 



Pi = 2 -2 ] J' J'. =0 for i = l,...,6, 

2 -2^^41 = 18 for i = 1,2,3. 

{ai,ai} 



Proceeding along the route of (|4.q ) reveals that for other Kovalevski exponents z = p{p—l) 
parameters satisfy the following equations: 

363^2-222-52/3 = 0, 2^-42-1/192 = 0. 

Both of them lead to complex Kovalevski exponents irrespectively of Pfad or concrete 
Bianchi type concerned. 



7 Conclusions 

Qualitative analysis of the dynamics of homogeneous Bianchi models revealed a remnant of 
the hierarchy of invariant manifolds. In the Bogoyavlenskii [p!^ ] formalism different Bianchi 
types lie on the components of the boundary of compact manifold on which respective 
dynamical system is defined. This boundary has its own boundaries thus generating a 
hierarchy in the class of homogeneous cosmological models which could be put into one- 
to-one correspondence with respective Bianchi type. The most general of homogeneous 
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models i.e. the so called mixmaster models are nonintegrable (e.g. [10, 11, |12[). We have 
therefore asked the question whether there exits a respective hierarchy of nonintegrability 
in the full class of homogeneous cosmological models. In this paper we answered this 
question at the level of algebraic integrability by showing that Bianchi VIq and VIIq 
models are algebraically nonintegrable since they have complex Kovalevski exponents. 
This property is propagated into more general Bianchi types (IX and VIII) according to 
the structure of invariant submanifolds. We have shown algebraic integrability of Bianchi 
type I and II models. 

We have also discussed the influence of matter terms in the energy-momentum tensor 
on the property of algebraic integrability of homogeneous cosmological models. It turned 
out that inclusion of dust, cosmological constant or radiation introduces algebraic nonin- 
tegrability in the system. One can wonder whether the cosmological constant for example 
should appear in the Kovalevski exponents and then give integer numbers in the limit as it 
tended to zero. However it is well known that one cannot obtain solutions of the Einstein 
equations without cosmological constant from the respective ones with cosmological con- 
stant by taking the limit A — > 0. This property is reflected also at the level of Kovalevski 
exponents. 

The literature concerning various aspects of homogeneous cosmological models is very 
big and diverse, e.g. Kramer et al. [|l^ gives a compendium of exact solutions of the 
Einstein's equations (and in large part is devoted to homogenous models), also a dynam- 
ical systems approach initiated by Bogoyavlenskii has considerably been developed by 
Wainwright and Ellis (for the review see 111]) and chaotic aspects of Mixmaster homoge- 



neous models have recently been investigated (e.g. |19|) in more detailed way. However, 
the existence of an exact solution does not imply the existence of algebraic first integrals. 
Similarly, the hierarchy revealed by the dynamical systems approach concerns the general- 
ity of solutions not the property of their integrability. Also the works revealing dynamical 
complexity of some models (in terms of fractal basin boundaries |l^) can by no means 
be understood as a rigorous proof of their nonintegrability. Therefore our results pre- 
sented above can be perceived as complementary to these existing in the literature and 
are another step toward deeper understanding of mathematical properties of homogeneous 
cosmological models. 
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